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A class of optimal quadrature formulas defined by V. L. N. Sarma in a probabilistic context is 
shown to be identical with a class of formulas defined previously, in a different manner, by P. J. Davis. 
As a result a contrast is drawn between the maximum of the error of a quadrature fomula over the unit 
sphere in a certain function space, and its average — in Sarma's sense — over that sphere. 
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A few years ago V. L. N. Sarma [1J l defined a probability measure on a certain space of func- 
tions and proposed that quadrature formulas be studied in terms of their average error on this 
space. If 

G(/)»G»(/)=2«r/(*r) (1) 

r= i 

is any /V-point quadrature formula, and E Q (f) =I(f) — (?(/) where /(/) is the integral which Q 
is approximating, we shall denote by "o-(£^)" the "RMS" average of E Q with respect to Sarma's 
measure. Sarma calls a formula Q\ "completely optimal" if <r(E Q \) ^ cr(E Q 'x) for every /V-point 
formula Q' N ; and he calls it "optimal in the weights" if cr(E Q <\ ) ^ a(E Q k) for every /V-point Q' N using 
the same points x\, x>?, . . ., x\. It is the purpose of this note to point out that, in the case of 
integration of functions of one variable, Sarma's optimal formulas are identical with those defined 
earlier by P. J. Davis in a quite different manner; and to discuss some consequences of this identity. 
In the one-variable case, Sarma's space consists of all those functions analytic in the open 
unit disk whose Taylor series 



satisfy the condition 



The integral considered is 



f(z) = £ a»z» (2) 



2l««l = i. 0) 



Kf) = j' i f(x)dx; (4) 
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and he showed that the average error a is given by the formula 

CT 2(£Q) = 2 3-'- 1 (EQ(x»)) 2 (5) 

an, p. 6io). 

In 1953, P. J. Davis [2] introduced another measure of quadrature error for analytic functions. 
In the notation of [3], Davis considers the integral (4) for /belonging to the Hilbert Space H 2 (Cr): 
the space of all those functions analytic in the disk \z\ <R (where R is any number > 1) whose 
Taylor series satisfy the condition 

J \a n \ 2 R 2n <™. (6) 

The inner product is defined by 

(f,g) R = 27rR Ja n W? 2n =[ f(zMz)ds (7) 

where a n and b„ are the Taylor coefficients of /and g respectively and C H is the circle \z\=R; 
s denotes arc length. For any Q of the form (1), with quadrature points x r inside Cr, E q turns out 
to be a bounded linear functional, and its norm is given [3] by 

\\E%=^ J /*-«» (£«(*»))«. (8) 

2itK «=0 
A knowledge of this norm affords us the error bound 

|/(/) -£>(/) I «I|£WII; 

and Davis defined optimality of quadrature formulas in the same manner as Sarma, with the 
norm \\E Q \\ in place of <r(E Q ). 

Comparison of (5) and (8) shows that 

«r(£ Q ) = (fj£) 1 Vlv5; W 

since the two measures of E Q differ by a constant factor, it follows that the quadrature formulas 
optimal in Sarma's sense are precisely those optimal in Davis' sense for the space H 1 (C ^ ). 

The H 1 spaces have many convenient analytic properties, and the norms U^H/* have been 
studied for various (?'s ([4], [5]). In [3] I showed that for Q = G N — the TV-point Gauss-Legendre 
quadrature formula (for a description of these formulas see, e.g., [6], chap. 8) — we have the bound 



\\E (; N\\ R ^AV2N+lp- N (10) 

where 

p = 2R 2 -l + 2R(R 2 -l) 1 ' 2 (11) 



/d 1/2 + o- 1/2 \ 1/2 
A = 32eV*ir-*p-« ( P P p _ m ) . (12) 
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If R = V3, then p = 5 + 2 V6, and it follows that ||£ r; v|| vs goes to zero about as fast as (0.1010) v . 
or faster. In a recent paper [7], Sarma and A. H. Stroud calculated v(E g n) for N=2, 3, 4, . . ., 
20, and observed that <j(E g n)I(t{E g n-\) was very close to 0.10. Because of the relation (9), this ratio 
is the same as 

\\e g »\U /ll£' ; -v-.|U, 

and these data seem to indicate that the left side of (10) does, indeed, give about the right rate of 
decrease for ||£ ,( ' ; a' , || /j . 

(Sarma and Stroud remark that (t(E g n)J(t(E g x-\) appears to be approaching 0.1 as N-^^. 
A closer look at their data suggests that the limit is somewhat higher. If one calculates the first 
differences of the successive values of cr(E ( 'x)lcr(E ( 'N-\) , it is seen that these differences are posi- 
tive and decreasing, and the ratios of successive differences are increasing. If this behavior con- 
tinues for higher /V, <t(E g n)I<t{E g n-\ ) must eventually become greater than 



(t{E G 20)I(t(E G 19)+A 20 

where 



tHtHth-- 



(13) 



A N =o-(E G N)l<r(E G N-i)-<r(E G N-i)/<T(E G N-2). 



The quantity (13) is greater than 0.1003; so it appears that the Sarma-Stroud data are consistent 
with the ratio of the cr's approaching (5 -(-2 V6) _1 ). 

It is interesting to compare Sarma's a with a very* directly related error bound. Following 
Sarma, let S be the space of all those functions analytic in the open unit disk whose Taylor series 
(2) satisfy 

J |fln|<°°, (14) 

n = 

and define 

\\f\\s = J \a„\. (15) 

w = 

S is a Banach Space, isomorphic to the space l\ of absolutely convergent series of complex numbers. 
Sarma's cr(E Q ) is a root-mean-square average of E Q over the unit sphere in S. In contrast, the norm 
of E Q — regarded as a functional on S — is the supreme of |£^| on the same unit sphere. This norm, 
which I shall denote by "||£^||s" allows the error bound 

|/Cfl-<?(/)l*P«IU-H/IU. /<£• (16) 

In the spirit of Davis' approach, ||£^||s would be used in place of a(E Q ). Let us compare the sizes 
of these two quantities: 

For every R > 1, if/e// 2 (Cfl) then/eS; and if we write 



||/||«=(/,/)j/ 2 = 2tt/? X l«n| 2 ^»| ,/2 , (17) 

^ n = 

we may compare \\f\\ R and \\f\\ s . Setting b„ = \a n \R n , we have 

2 \a n \) = 2M- 



S 6 *- ){Z R - 2n 

n = o / \ n = 

87 



R 1 



R 2 

X «_ A 



by the Schwartz inequality; and it follows that 

11/11,. =s KiRMU 
where 

K(R) = ' 
Therefore 



R Y' 2 



2n(/? 2 -i) 

\E%f) 



\n = sup 

AS 11/ ||S 



(18) 
(19) 






\EHf)\ 

K(R)feH 2 (C r ) ||/||« 



sup 



K(R) 



\\E%- 



(20) 



The experimental evidence seems to indicate that, for Q = Gy, the asymptotic behavior of ||^ 6 ^||ii 
for large TV is correctly given (at least to within factors like N a ) by the right side of (10). If this is 
so, then (20) shows that \\E (; 4 S # 0(C~ N ) for any C> 1, while (by (9), (10), and (11)) a(E (: x) 
= 0(3 _A 7V 1/2 ). Equation (16) may be a very conservative bound; in numerical experiments [8] an 
analogous bound has been found to be considerably larger than the actual quadrature error. Since 
Sarma's "average error" cr(E (; x) appears to be very much smaller than the right side of (16), it 
would be interesting to know whether it does, indeed, represent the "average" behavior of quad- 
rature formulas for the integrands arising in the course of practical numerical work. I do not see 
how to answer this question. 



References 



[1] Sarma, V. L. N., Eberlein measure and mechanical quadrature 
formulae. I: Basic theory. Math. Comp. 22, 607-616 (1968). 

[2] Davis, P. J., Errors of numerical approximation for analytic 
functions. J. Rat. Mech. and Anal. 2, 303-313 (1953). 

[3] Haber, S., The error in numerical integration of analytic func- 
tions, to appear. 

[4] Hammerlin, G., Uber ableitungsfreie Schranken fiir Quadratur- 
fehler,Num. Math. 5,226-233 (1963). 

[5] Barnhill, R. E., Asymptotic properties of minimum norm and 
optimal quadratures, Num. Math. 12, 384-393 (1968) . 



[6] Hildebrand, F. B., Introduction to Numerical Analysis (McGraw- 
Hill N.T., 1956). 

[7] Sarma, V. L. N., and Stroud, A. H., Eberlein measure and 
mechanical quadrature formulae. II: Numerical results, 
Math. Comp. 23,781-784 (1969). 

[8] Barnhill, R. E., Optimal quadratures in L 2 (£,>). I. SIAM J. 
Numer. Anal. 4,390-397 (1967). 



(Paper 75B1&2-348) 



88 



